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EFFECT OF EXTERNAL DISTURBANCES AND DATA RATE ON 
THE RESPONSE OF AN AUTOMATIC LANDING SYSTEM 
CAPABLE OF CURVED TRAJECTORIES 

Windsor L. Sherman 
Langley Research Center 

SUMMARY 

The automatic landing system described in NASA TN D-7611 is capable of guiding a 
large transport airplane on curved decelerating trajectories to a landing on the airport 
runway. The system which takes over when the airplane is close in, 6000 m or less from 
the landing point, will, for the acquisition of guidance data, interface with the microwave 
landing system. A nonstatistical study has been made of the effect of wind, shear, turbu- 
lence, data sample, and control-actuator natural frequency on the response of the system. 
The results indicate that the system functions well in the presence of wind shears and tur- 
bulence. However, for steady wind the guidance laws, particularly the turn algorithm, 
must contain proper compensation for the effect of wind. The system had satisfactory 
response for data sample rates down to five samples per second and control-actuator 
natural frequencies of 5 Hz. 

In limiting cases, those cases where the roll angle reaches zero as the airplane 
touches down, the gains in the turn algorithm computed by the methods given in NASA 
TN D-7611 produced unsatisfactory landings when winds and turbulence were included. A 
revised method for determining the turn -algorithm gains, presented herein, computes turn- 
algorithm gains that give satisfactory landings in the limiting cases. The use of gains 
determined by the methods in nonlimiting cases also improved the system response. 

INTRODUCTION 

Reference 1 presents the results of a study of a possible automatic landing system 
in a bland environment. The system studied is capable of guiding large jet transport -type 
aircraft over steep, curved decelerating trajectories over the portion of the landing oper- 
ation that occurs just before touchdown. It takes over when the airplane is 6000 m or less 
from the landing point. The data requirements of the automatic landing system are com- 
patible with the data output of the microwave landing system (MLS). This report extends 
the results presented in reference 1 to a nonbland environment. The effects of steady 



wind, wind shear, random turbulence, data sample rate, and control-actuator natural fre- 
quency on the response of the airplane and system are shown. In the case of wind shear, 
the heading and/ or the magnitude of the shear was dependent on altitude. In addition to 
the overall response study, the effect of initial conditions, position, altitude, and environ- 
mental conditions on the gains in the turn algorithm was studied. The present study, which 
is analytical in nature, was made on a large digital computer for the airborne automatic 
landing systems described in reference 1, 

Ideally, the description of automatic landing system performance in a nonbland 
environment would be a statistical study involving multiple landings at various initial con- 
ditions so that rms errors could be determined. However, such a study involves rather 
complete knowledge of the system hardware characteristics as well as the specification 
of the external disturbances. In addition, a study of this type consumes vast amounts of 
computer time. Inasmuch as the hardware for the autoland system described in NASA 
TN D-7611 is, not known and in order to conserve computer time, typical initial conditions 
including those at the limit of the system performances were studied on a nonstatistical 
basis. The results obtained, while not suitable for the determination of operational limi- 
tations, demonstrate the response of the system for various data sample rates and atmos- 
pheric disturbances, and show the effect of the initial condition on turn -algorithm gains 
together with a method for determining these gains. 

SYMBOLS 

The International System of Units (SI) is used throughout this paper. All angles 
are measured in radians. The coordinate systems referred to are discussed in 
appendix A. 

Baa radar azimuth angle in R^p coordinates 

Bag desired radar azimuth angle in Rjp coordinates 

C, rolling- moment coefficient due to sideslip 

H altitude 

K desired angle between runway centerline and velocity vector at point B of 

figure 9 

k^ variable gain for the pseudo-radar azimuth angle error of the turn-control 

algorithm 
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variable gain for the total error used to determine roll command 


= 0.3 

gain used in computation of k.2, may be variable 
pitch rate 

projection of range vector in the z^,Z2 plane 
special radar coordinates used in calculation of 
radar range components in Zj coordinates 
component of ground speed in direction 
windspeed along Rjp axis evaluated at point A only 
windspeed along rj^-axis 
airplane speed 

windspeed evaluated at point B only 
windspeed along 
windspeed along rjg-axis 
vertical speed in Zj coordinates 
radar coordinates 

airplane coordinates referred to principal body axes 
quasi-inertial coordinates 
angle of attack 
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|3 

4 


y 

a. 




Vi 


^AA 

®AC 

®AW 

®AAV 

®ACV 

= 0 
^ c 


sideslip angle 

sideslip angular rate, rad/sec 
flight -path angle 
aileron deflection 
rudder deflection 

runway coordinates, referred to the desired landing point 
inertial coordinates, translate with airplane 
pitch angle 

radar azimuth angle in Zj coordinate 
desired radar azimuth angle in Zj coordinate 
wind correction angle 

radar azimuth angle in z^ coordinates referred to velocity vector 
desired radar azimuth angle in Zj coordinates referred to velocity vector 
roll angle 

heading angle, airplane centerline 
heading angle of runway 
heading angle of velocity vector 
special heading angles used to determine 
dummy yaw angle used in determination of kj 
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Subscripts : 


c command 

i initial 

j coordinate index, j = 1, 2, 3 

A dot over a variable indicates differentiation with respect to time, 

THE AUTOMATIC LANDING SYSTEM 

The automatic landing system used in this study was fully described in reference 1. 
Briefly, data from an inertial table and either an airborne radar or the microwave landing 
system are processed in airborne computers to obtain guidance signals that direct the air- 
plane along steep, curved approaches to the landing point. A moderate speed change, about 
5 m/sec, was accomplished during the landing maneuver. The system operates diming final 
approach and landing, assuming control of the airplane between 3000 m and 6000 m from 
the landing point. No major changes were made in the system as a result of the introduc- 
tion of cross winds, turbulence, and sample rate. However, some modification of the head- 
ing angle and radar -azimuth -angle feedback were required to handle cross winds. In the 
presence of head winds, a change to the flare computer also was needed in order to insure 
that the airplane had the proper attitude angle at touchdown. These changes will be dis- 
cussed when the appropriate results are presented. The block diagram of the autoland 
system with the changes found necessary for atmospheric disturbance is presented in 
figure 1. 

The airplane that was modeled for this study was a large four -engined jet transport 
that had a mass of 90 719 kg. Complete data for the airplane are given in reference 1. 

Disturbances Used 

Steady winds with and without shear and turbulence were introduced as atmospheric 
(external) disturbances. In the case of steady wind, the wind was always considered to be 
parallel to the ground plane. Shear was introduced on the steady wind by two methods: 

(1) a magnitude change with altitude and (2) a direction change with altitude, which is 
referred to rotation. The magnitude changed at the rate of -0.048 m/sec per meter, from 
a windspeed of -25.80 m/sec at a 540-m altitude, and the rotation was 0.0058 rad/m. 

The turbulence used in this study was based on a random-number program, the out- 
put of which is a normal distribution of random numbers between ±1.0. As long as the 
initialization is not changed, the same set of random numbers is generated. The random 
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numbers generated by the program were used to generate an acceleration which was 
integrated to give a speed. The turbulence was generated in inertial coordinates, the 
Tjj-system, and transformed into airplane body axes. A separate random-number gener- 
ator with different initialization was used for each of the Tj^-axes. Three levels of turbu- 
lence, mild, medium, and strong, were generated. The rms gust speed and the maximum 
gust speed for these three levels of turbulence are given in table I. Because the initializa- 
tion of the random-number program was not changed during the course of the investigation, 
the turbulence repeated from run to run. Typical time histories of the medium turbulence 
are shown in figure 2. As can be seen from the time history, the turbulence consists of 
continuous gusts of relatively high intensity whereas atmospheric turbulence generally 
consists of brief periods of large gusts separated by periods of relative calm, A so-called 
patchy turbulence was obtained by using the output of one of the random -number generators 
to control the intensity of the gusts along the three ■rj^-axes. This produced a turbulence 
(see fig. 2) time history (labeled patchy) that more closely resembles atmospheric condi- 
tions than does the unmodified medium turbulence. The rms gust speed and the maximum 
gust speed are given in table I, During the study of the effect of turbulence on the system 
response all of the turbulences listed in table I were used. However, when combined with 
other effects, such as steady wind, either the medium or patchy types of turbulence were 
used. 

In addition to the external disturbances, the sample rate ^the number of times per 
sec that the continuous output of the guidance computers was read to obtain new values of 
<p and q ) and the control-actuator natural frequency were varied. Five sample rates 
between 1000 samples per second and 5 samples per second were used. The control actu- 
ators programed for the autoland system were considered typical of a modern commercial 
jet airplane and had a natural frequency of 30 Hz. Natural frequencies as low as 5 Hz were 
used to determine the effect of varying this parameter on the system response. 

Acceptable Touchdown Conditions 

In the investigation reported herein it was assumed that the landing was taking place 
on a runway 3000 m long and 50 m wide. The desired touchdown point was on the center- 
line 100 m from the end of the runway. The values of <p, !//, ?2? ^.nd y, because 

of the reference used, indicate the errors from the ideal touchdown condition in which all 
these variables would have a value of zero. In this study, for the purpose of determining 
if a landing was satisfactory, an arbitrary set of conditions was established. If a landing 
fitted within the following limits, it was said to be satisfactory: 

cp = ±0.06 rad 

\jy = ±0,01 rad (except for steady winds) 
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= +500 m 

^2 = ±10 m (based on airplane model used in study) 
0 > 7 ^ -1.98 X 10"^ rad 


0 < w„. = 1.0 m/sec 
0.017 ^ e ^ 0.061 rad 

The heading angle reported in the tables is the heading of the airplane centerline with 
respect to the runway. The autoland system controls the pointing of the velocity vector, 
so in the presence of wind the reported heading angle is that which is required to direct 
the resultant ground-speed vector down the runway centerline. This approach to control 
minimizes speed perpendicular to the runway centerline. As the heading angle is a func- 
tion of windspeed and direction, no real tolerance criterion is possible. However, the 
values of given in table II for cases with wind place the velocity along the centerline 
within reasonable limits. 


RESULTS AND DISCUSSION 

The results obtained in this study can be divided roughly into the effect of wind con- 
dition and the effect of system conditions. Table II summarizes some of the results 
obtained in this study. Data are given for two important points in the landing maneuver, 
the start of the flare and the actual touchdown. The nominal point at which the flare is 
started is an altitude of 20 m. For the data given in table II and all tables that follow, the 
start of flare occurs in a narrowband about this nominal altitude, usually within ±10 per- 
cent, Case 1 in table II is the data for the autoland system described in reference 1. 
There were no disturbances in the model when these data were obtained. Because the 
results presented in reference 1 showed that if the gain was set correctly, the 
response of the system was not initial -condition dependent, only one initial condition 
was studied in detail, the one given in table II. 

Effect of Turbulence 

Cases 2 to 4 of table II show the effect of turbulence on the autoland system. An 
inspection of the ground and altitude tracks showed that there v/as little effect on them 
from the turbulence. The touchdown conditions with the exception of 0 and Wzj were, 
in general, acceptable. 
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The principal effect of turbulence was on the motion of the airplane itself. The 
effect on the airplane can be seen by comparing time histories of aircraft motion during a 
typical landing maneuver presented in figure 3, None of the motions presented appear too 
drastic and, as can be seen, the patchy turbulence produces slightly less effect than the 
medium turbulence. It should be emphasized that the autoland system used in the turbu- 
lence studies was the same system described in reference 1. No gains or other changes 
were made to the system. In general, the turbulence had no important effect on the ability 
of the system to execute a landing. During the remainder of the study when turbulence was 
used, it was the medium or patchy turbulence. The light turbulence had no effect on the 
-system. The autoland system was able to control the airplane in the heavy turbulence; this 
was considered an extreme case because of the high gust speeds and high touchdown speed. 


Steady Winds 


Steady winds with speeds of -25.80 m/sec and 14.007 m/sec were used at various 
angles to the runway centerline between 0 and 7t/2. The ground tracks for the basic sys- 
tem (case 1) and the basic system with wind (cases 7 and 9) are shown in figure 4. The 
windspeed was -25.80 m/sec at tt/ 4 rad to the runway centerline. With the wind present 
the ground track deteriorated as shown in figure 4 and touchdown conditions were not 
acceptable (see case 7, table II). The altitude track was also unacceptable. This deteri- 
oration was due to the fact that the guidance system basically controls the direction of the 
inertial velocity vector but the controlling parameters ■ and 0^^ are referenced to 
the airplane centerline. 


The linear velocities and along the R 2 ^-axes were used to determine a 

wind correction angle ®AW which is given by the following expression: 


^AW 


= tan 



( 1 ) 


The angle 0^-^ was combined with the heading angle of the airplane to obtain the 
heading angle of the velocity vector ^ + ®Aw)- In addition, through the use of this 

angle, Rz^ and Rzg were transformed so that Rzj^ was along the velocity vector and 
Rz 2 perpendicular to it. The transformed values of Rzj and Rz 2 which are R^^ 
and R 2 p were used to compute two new radar angles 0^Qy and 0^^y that are refer- 
enced to the velocity vector. These new angles were used in the turn algorithm to deter- 
mine cp . The modified equation for cp is 

U C/ 




4/ - U/ - 


^ 3 ^ + ki(0^Cv " ®AAv) 


( 2 ) 


This change of reference from the centerline to the velocity vector restores the accuracy 
of the autoland system (see case 9 in fig. 4). Cases 8 to 11 of table II summarize the 
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results obtained for steady winds and a comparison of the ground tracks is presented in 
figure 5. No changes in the letdown guidance were required and there were no effects on 
the altitude tracks. 

The changes to the turn algorithm, summarized in equation (2), introduce wind and 
turbulence information into the determination of the turn commands. Cases 5 and 6 of 
table II were run in order to determine if turbulence affected the modified turn-control 
system. The data presented in the table for these cases indicate no effect on the system 
response. There were no significant changes in the ground and altitude tracks for these 
cases. 

Two additional control changes were introduced to improve control with steady wind. 
These were and 6^ feedbacks to the ailerons and a change in the flare computer. 

The former, which coordinates the turns, will be discussed first. Under some conditions 
of sideslip, the rolling produced by is sufficient to cause the airplane to roll in the 

wrong direction. By feeding back /3, and dj- to the aileron the undesired rolling 

moments were canceled. The commands from the guidance then modify the corrective 
6„ to obtain the desired roll angle. Because the correction is dependent on /3 the cor- 
rection is self-canceling as it approaches zero and as /3 approaches zero. This type of 
correction is not considered a mandatory part of the autoland system as the necessity for 
it depends on the aerodynamics of the particular airplane. The details of the change are 
shown in the block diagram in figure 1. The second change was included when an inspec- 
tion of the touchdown conditions showed that the flare routine was not pitching the airplane 
sufficiently to insure a safe landing. A nonstandard procedure was adopted to obtain the 
desired pitch angle. The flap controls were modified so that the flaps were retracted until 
the airplane reached a pitch angle of 0.052 rad, a pitch angle that gives a reasonable pitch 
attitude at touchdown. The retraction was started at the time the flare was started. 

In addition to steady winds, vertical wind shears and winds that changed directions 
were used. The introduction of these wind conditions did not cause the system to abort 
landings (see cases 12 to 14 of table II)’, nor did the airplane and system response (i.e., 
ground and altitude tracks) deteriorate. 

When the wind is not directly along the runway, the airplane is flown so that the 
resulting velocity vector points in the correct direction; thus as the airplane approaches 
the runway, it is yawed, or crabbed, with respect to the runway centerline. Much has been 
done on the decrab maneuvers, that is, the alinement of the airplane and runway centerlines 
just before touchdown. Reference 2 reports a typical decrab -maneuver study for transport- 
type aircraft. The results given in reference 2 indicate that a decrab maneuver is feasible 
and that the systems used could be incorporated in the automatic landing system of this 
study. For these reasons no extensive study of the decrab maneuver was made; however. 
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enough work was done to establish that conditions at flare initiation and touchdown would 
make it feasible to incorporate a decrab maneuver. 

The next logical step was to combine the steady winds and the turbulence. Ground 
and altitude tracks are shown in figure 6 for medium turbulence and a steady wind of 
-25.80 m/sec; this is case 15 of table n. As can be seen by a comparison of figures 5 
and 7 the combined wind and turbulence caused little change in the tracks. The airplane 
responses were changed very little from those given in figure 3 for the medium turbulence. 
Both the medium and patchy turbulence were used in these runs; neither caused unsatis- 
factory results. Case 16 of table II combines all types of wind effect that have been used 
in this study. As can be seen, this combination did not have an adverse effect on either 
the touchdown conditions or the gromid and altitude tracks. 

Sample Rate 

In the sample-rate studies it was assumed that the outputs of the angle and distance 
sensors were continuous and that the guidance computers used a finite number of samples 
of data per second for determining the guidance commands. Sample rates of 1000, 100, 

50, 10, and 5 samples per second were used. The data sample rate was also the refresh 
rate for guidance commands. If 10 samples of data were taken per second, the guidance 
computer output per second was 10 discrete commands; one for each sample of data. The 
sample-rate tests were run with medium turbulence and are recorded as cases 5 and 17 
to 20 of table II; case 5 is the basic case. 

Data for these cases, presented in table II, show that there was little change in lat- 
eral variables with sample at the start of flare and the touchdown conditions were within 
acceptable limits. Ground tracks for case 5 (1000 samples per second) and case 20 
(5 samples per second) are shown in figure 7. As can be seen, there are only small dif- 
ferences in the ground tracks of these extreme cases. Inspection of the results of the 
vertical mode indicated that up to the initiation of the flare the results were similar to 
those for the horizontal control. Conditions at the start of flare indicated that an accept- 
able flare could be performed. However, at touchdown the only variable that at all times 
was within acceptable limits was The vertical velocity at touchdown increased from 

4.48 X 10" 1 m/sec for 1000 samples per second to 2.30 m/sec for 5 samples per second, 
which is unacceptable. Oddly enough the pitch angle at touchdown changed from an unac- 
ceptable value of -6,46 x 10 at 1000 samples per second, to an acceptable value of 

O 

3.92 X 10”'^ at 5 samples per second. The increase in the vertical touchdown speed Wzj 
is accounted for by a failure of y to decrease toward zero, as the sample rate decreased. 
Because | y | increased with decreasing sample rate, a refresh rate between 50 and 
100 commands per second is indicated. The most practical way to accomplish this is to 
add a predictor to the flare computer so that commands between data samples can be 
generated. 
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Case 21 of table II combines the steady wind of -25.80 m/sec with the medium turbu- 
lence at 10 samples per second. The flare conditions for this indicate a good touchdown 
can be achieved but, as with the other 1000-samples -rate cases, 9 and Wjjj proved to 
be imacceptable. 

The fact that no uncorrectable deterioration of the system occurred at 5 and 10 sam- 
ples per second is most important, as these are the data rates for the MLS system in 
category II and III conditions (see ref. 3). Because of their relationship to the MLS, sample 
rates of 10 and 5 samples per second were used for further study. No increase in the com- 
mand refresh rate was incorporated in the system. 

Control-Actuator Natural Frequency 

The control actuators modeled for the autoland system were representative of the 
control surface actuators of large modern transport aircraft. Because this type of actu- 
ator has a natural frequency of about 30 Hz and has damping ratios greater than one, a 
first-order representation for the actuator was used in the study. In order to gain some 
idea of the effect of actuator natural frequency on the system response, a run was made 
with actuators with a natural frequency of 5 Hz. The sample rate was 5 samples per sec- 
ond. The ground track for this case differed very little from that for case 20 (see fig. 7), 
and the altitude track was satisfactory to the start of the landing flare. The touchdown 
conditions are given in table H as case 22, and are satisfactory except for 9 and Wg.. 
These results were typical when the natural frequency of the three control-surface actu- 
ators was the same. However, if the natural frequency of each actuator was different 
the system response was much poorer than that for case 22, The worst response, which 
was not acceptable, occurred when the aileron and rudder actuators had the same natural 
frequency, and the natural frequency of the elevator actuator was lower than the others. 

Discussion of Results Presented in Table H 

In cases 1 to 22, when properly compensated for environmental conditions, the lateral 
(turn control) guidance was able to reduce the lateral displacement errors and heading- 
angle errors to acceptable values. The absolute value of the roll angle at touchdown was 
equal to or less than 5.35 x 10“^ rad which is within the acceptable limits. This was not 
the case in longitudinal guidance, as for 13 of 22 cases in table II either 9 and/or Wgjj 
did not fall within acceptable limits. Of these 13 cases, there were 3 (cases 3, 19, and 20) 
where the value of only Wg^ was unacceptable. In cases 6, 12, 14, 15, 16, 17, and 18 all 
of the longitudinal touchdown conditions were satisfactory. The average values of three 
longitudinal parameters, the pitch angle 9 , the flight -path angle y, and the sink rate Wzj 
were computed at the start of the flare and are 
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9 = -6.83 X 10"^ rad 
y = -5.48 X 10"^ rad 
W 2 J = 4.023 m/sec 

In all cases where one or more of the longitudinal parameters did not meet the criterion 
for acceptable touchdown conditions, one or more of the parameters 6 , y, and Wz^ 
deviated from the above values at flare initiation. For instance, in case 1 where only W 2 j 
did not meet acceptable touchdown conditions, the sink rate at the start of the flare was 
1.07 m/sec above the previously given average. This suggests that more precise control 
of the state conditions at the start of the flare would help to obtain acceptable touchdowns. 
One interesting point connected with the acceptable landings is that there is either no 
steady wind in the problem or a wind with a shear pattern that reduces the windspeed to 
zero at ground level. This implies that there is an effect of steady wind on system per- 
formance that has not been completely identified. 

Limiting Initial Conditions 

In reference 1 limiting initial conditions were defined as those cases in which the 
roll angle reached an acceptable value as the airplane touched down on the runway. One 
such case has the following initial conditions: 

= -3000.0 m 
?2 = -4000.0 m 
H = 540 m 
= 0.0 

When this initial condition was used in reference 1 a satisfactory landing resulted. The 
results for the present case with a steady wind of -25.80 m/sec at tt/4 rad to the runway 
centerline and patchy turbulence are given as case 23 of table in. All touchdown conditions 
were unacceptable and the conditions for 9 , y, and W 2 j at flare initiation differed greatly 
from the previously given average values at this point. The complete unacceptability of 
this case is borne out by the ground track (see fig. 8). 

A rerun of this case without the steady winds showed large improvement in the touch- 
down conditions. This result again pointed to an effect on the system of a steady wind. 
Although the determination of the parameters in the lateral guidance had been changed to 
account for steady wind, the gains k^ and kg were still being computed by the methods 
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given in reference 1. As was pointed out in reference 1, system performance is very 
sensitive to these gains and these gains are sensitive to initial conditions; therefore, the 
logical place to look for an uncorrected wind sensitivity was these gains. 

The Gains and ^2 

The gains k^ and k 2 are used in the lateral-guidance turn algorithm (eq. (2)). 

The gain k^ is recomputed whenever the direction of the turn is changed and the gain kg 
is continuously calculated as a function of the heading angle. For this method k^^ calcu- 
lated at the beginning of the landing maneuver does not take atmospheric conditions into 
account; however, any subsequent calculation of kj does take these conditions into 
accbunt. Therefore, it is the first calculation of k^ for which a revision was made. 

The new method of calculating the initial value of k^ was changed to consider winds. 

(See appendix B.) A comparison of the results for kj^ computed by the method of refer- 
ence 1 (case 23) and that of appendix B when equation (B9) is used to calculate kj 
(case 25) is given in table III. Cases 24 and 26 show the effect of small changes, approxi- 
mately 2.5 percent in the kj^ of case 25. The change of kj to 1.759 produced a set of 
satisfactory landing conditions. While the values of 0, y, and Wzj at flare initiation 
are not close to the previously given average values of these parameters for the acceptable 
cases of table II, they lie within the spread of the data of these cases. Increasing k^ by 
approximately 2.5 percent produces a completely unacceptable landing, case 24. Decreas- 
ing k^ by about the same amount improved the values of the parameters cp, \j/, and 
but the overshoot in the direction, about 800 m, may not be acceptable for short run- 
ways. Case 25 was adopted as a base case for the study of the effect of the gain kg. Fig- 
ure 8 compares the ground tracks for cases 23 and 25. 


The Effect of Varying the Gain kg 

The gain kg appears in equation (2) as a multiplicative factor of the total error and 
converts this heading angle and position error to a bank-angle command. The gain kg is 
given by 


k<j = — e 


4 .-I'^c-'^vl 


( 3 ) 


where k 4 is a constant. The only way to vary ^2 is to change k 4 . These changes 
in k 4 are listed in table IV. However, because k 2 appears in equation (2), the changes 
are referred to as variations in k 2 . The changes in k 2 do not affect kj so no values 
of kj are given. The values of kj are the same as corresponding cases of table III. 
Two basic cases were used for table IV: case 23 in which the landing was unsatisfactory, 
and case 25 in which the landing was satisfactory. An inspection of cases 23, 23a, 23b, 
and 23c shows that when k 2 was not correct, variations in kg did not improve the 
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landing conditions but actually shows further deterioration in some of the parameters 
such as the roll ai^le cp. However, when kj had a value that gave a satisfactory 
landing (case 25) changes in k 2 definitely affected the landing. When k 4 was 
decreased from 4.8 (case 25, table IV) to 4.4 (case 25a, table IV) all of the landir^ 
parameters improved, cp and ^2 were reduced by about an order of magnitude with 
smaller improvements in the overshoot in the direction, and there was a small 
improvement in the pitch attitude at touchdown. When k 4 was reduced still further 
to 4.0, the landiiig conditions moved back closer to those obtained when k 4 was equal 
to 4.8. Accordingly, a value of 4.4 was adopted for k 4 . The values of 0, y, and Wzj 
at flare initiation did not match the average values of these parameters for the acceptable 
cases of table II and the only one that fell within the spread of the data was Wzj. These 
results appear to indicate that correct values of kj and k 2 are the important factors 
in achieving a successful landing. 

To this point the only value of used in the study of the gains is 0.0. Other 
values of were tried and it was found that a small correction to the gains kj and 
k 2 based on initial heading angle was required to obtain the best landing conditions over 
a wide range of 4^^, The final equations used to compute k^ and k 2 are equations (BIO) 
and (B14) of appendix B. 

Table V summarizes the conditions at the start of flare and at touchdown and ground 
tracks are shown in figure 9 when equations (BIO) and (B14) were used to calculate kj 
and k 2 < These results are for the following initial conditions: 

-6000 m, ? 2 = -4000 m 
H = 540 m, 4/^ = 0.0, tt/4, n/2 
= -3000 m, ?2 = -4000 m 
H = 540 m, 4^^ = 0.0, tt/4, tt/2 

with patchy turbulence and a steady wind of -25.80 m/sec at tt/4 rad to the runway center- 
line. The sample frequency was 10 per second. The gains k^ and k 2 which appear in 
the turn algorithm primarily affect the ground track of the airplane; however, because the 
flight-path command y^ is determined in a plane that is perpendicular to the ground 
plane and passes through the start-of-flare point and the instantaneous position of the air- 
plane, some effect of the lateral guidance can be expected on the letdown of the airplane. 

As can be seen from table V and figure 9, the lateral control is very good and the airplane 
is well alined with the runway before touchdown with very small roll angles and lateral 
displacements from the runway center. The body yaw angles indicate that the pointing 
of the velocity vector with respect to the runway centerline is also satisfactory. The 
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longitudinal touchdown conditions were not satisfactory with regard to pitch angle 
(cases 27a and 27b) which was too small; in case 27b, W 2 j is slightly high. Case 28c 
is also interesting. The pitch angle and Wzj are both good although 0 is a little high, 
the large discrepancy being the overshoot of the touchdown point. In this case, the over- 
shoot is over 800 m but is still within the first third of the assumed 3000-m runway which 
is acceptable. A detailed study of the results indicated that the control system was func- 
tioning in the intended manner, that is, the flight-path angle was approaching zero and the 
airplane was pitching up to the proper touchdown attitude. This lack of proper perform- 
ance in the flare region is in marked contrast to the performance of the longitudinal control 
-system in the letdown prior to the flare. In this preflare region the letdown was well con- 
trolled with end conditions that should have permitted the execution of an acceptable flare. 

The inconsistency of flare performance occurred not only in the results presented in 
this section, but in the results presented in the other sections of the paper. In addition to 
command refresh rate during the flare, analysis of the data indicates that the slowness of 
the speed control and/ or the lack of coordination between the flight-path control system 
and speed control system may be the causes of the flare -touchdown problems. A better 
coordinated flare controller or a decoupled control system for the longitudinal mode which 
controls the flight path, pitch angle, and speed might help these problems. Research on a 
decoupled control system, steady -state decoupling only, for STOL aircraft (refs. 4 and 5) 
has shown this type of system to be highly effective for precise control of the longitudinal 
mode. 


CONCLUDING REMARKS 

An automatic landing system capable of guiding an airplane over steep, curved paths 
to a landing has been studied to determine the effect of wind, turbulence, data sample rate, 
and control-actuator natural frequency on the ability of the system to execute a landing. 

It was found that wind, including shears, and turbulence did not impair the ability of the 
autoland system to execute a landing. The gains in the turn algorithm, kj (variable gain 
for the pseudo-radar azimuth angle error of the turn-control algorithm) and kg (variable 
gain for the total error used to determine roll command), were critical factors in a suc- 
cessful landing and should be determined by the method presented in this paper. There 
were no observed effects of changing the control-actuator natural frequency as long as all 
three control actuators had the same natural frequency. Sample rates down to five samples 
per second were studied. At microwave -landing-system sample rates of 5 and 10 samples 
per second good system response was maintained and successful landings were executed. 

The lateral response, alinement, and centering of the airplane with the runway were 
consistent and precise, as was the longitudinal response to the start of the flare. After 
the start of the flare there was a lack of consistency in the longitudinal response resulting 
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in unacceptable or marginal values of touchdown distance, pitch angle, flight -path angle, 
and vertical speed at touchdoivn. The most probable causes were a low command refresh 
rate in the flare, the slow response of the speed control, and a lack of coordination between 
the speed system and the flight -path control system. A steady-state decoupled control sys- 
tem, that controls speed as well as flight path, should be investigated as a longitudinal con- 
trol system. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
May 30, 1975 
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APPENDIX A 


COORDINATE SYSTEMS USED IN STUDY 
Five coordinate systems were used in the autoland system: 


(1) inertial coordinates 


(2) quasi-inertial coordinates 

^3 

(3) airplane coordinates 

^3 

(4) radar coordinates 

^3 

(5) runway coordinates 

^3 


Figure 10 shows these coordinate systems and their relationship for a typical landing 
situation. The subscript j in the symbol for the coordinates takes on the values 1, 2, or 
3, the values denoting a specific axis. When j = 1, the axis corresponds to the x-axis of a 
usual x,y,z coordinate system. Similarly, j = 2 corresponds to the y-axis, and j = 3 
the z-axis. 

The inertial coordinates do not rotate and are forced to translate with the airplane. 
The ?]j-axis has the same direction as the runway centerline. The quasi-inertial coordi- 
nates are rotated by the angle i//, obtained by a positive rotation about the rjg-axis. The 
z^-axis now points in the direction of flight. The airplane coordinates are related to the 
quasi-inertial system through the two additional angles 9, obtained by a position rotation 
about the Z2-axis, and xj/, obtained by a positive rotation about the y^-axis. 

The runway coordinates use the-touchdown point as a fixed origin and the is 
positive toward the far end of the runway. These coordinates, except for j = 3, are the 
negative of the inertial system. This system is used only for the reporting of the results. 
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APPENDIX B 


DISCUSSION OF THE GAINS kj AND kg THAT ARE USED 
IN THE TURN ALGORITHM 

The basic equation in the turn algorithm (eq. (2)) determines the bank-angle com- 
mand. Of the three gains in this equation k^^ and kg are the most important. These 
gains, through their influence on bank angles, influence the shape of the ground track and 
the letdown commands. 


The Gain k^ 

Figure 11 represents the ground track for an initial condition such as case 23 of 
table III. The points marked A and B in this figure are the points at which kj^ is cal- 
culated. The gain kj is calculated at point A so that the airplane will turn onto a head- 
ing that is approximately tt/2 rad with respect to the runway and have zero bank ai^le 
when at point B, The second value of kj is calculated at point B to give good landing 
conditions at the touchdown point and to prevent the airplane from making a reverse turn 
at point B. Experience has shown that these two conditions are compatible and the condi- 
tion applied is = 0.0 at point B. The basic modification required to the method given 
in reference 1 for the computation of kj at point A is the proper inclusion of the wind 
done in the following manner: 


^Ip 


cos (K - 

'P) + Rzg sin (K - \{y) 

^2p 

- -Hz 

j sin (K 

- x}/) + Rz2 - xl/) 


= u 

cos K + 

V sin K 

wp 

w 

w 


= -u 

sin K + 

V cos K 

wp 

w 

w 


(Bl) 

(B2) 

(B3) 

(B4) 


where K is a constant that is the angle between the runway centerline and the desired 

heading of the airplane at point B and u,„ and v,„ are the wind components along and 

w w 

perpendicular, respectively, to the runway centerline. The calculated parameters 
^2p’ ^wp’ ^wp to determine \}/^, a.nd B^^ as follows: 


= K + 


(B6) 
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^ = tan“^ \ 7^\ - i>/ 


AA 


IP/ 


Bag = tan 


-1 


R 


100.0 

Ip " Baa 


(B7) 


(B8) 


where 

eters 


Uq is the component of the ground speed in the 
Bag? ®AA’ '^5 were used to determine 


kf = 




®AC “ ®AA 


direction. The param- 
from the following equation: 

(B9) 


This value of kj proved adequate for an initial condition such as case 23 or case 1. 
However, when the initial heading angle is between 0.0 and tt/2, a modification based on 
heading angle was required. The final equation for kj which is adequate for all is 



^c -^5 
®AC “ ®AA 


- 0.09 sin(2xl/^ 


(BIO) 


The value of k^^ determined by equation (BIO) was used until the airplane reached 
point B in figure 11, at which point kj^ was recalculated. Point B was identified by the 
passage of the roll angle through zero, the actual condition being 


if (-0,002 ^ <p ^ 0.002) recalculate kj 

^ - ^3'^ 

1 ®ACV " ®AAV 


(BID 


The value of k^ determined by equation (Bll) was used from point B until the airplane 
touched down. 

The parameter K gives the desired heading angle between the runway centerline 
and the flight path at the point where k^^ is recalculated. The value assigned to K is 
Jx/2 and was used because it eliminates the possibility of singular conditions occurring in 
equations (BIO) and (Bll). 

The program used for the computation of kj contained a lower limit on the value. 
That is, if the calculated k^ was less than 1.483 for a left -turn approach or less than 
1.524 for a right -turn approach, it was set to 1.483 or 1.524. Case 28a of table V was the 
only case of this study where k^ passed below the limits. A study of case 28a showed 
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that better landings were obtained when the lower limit on the was dropped. Accord- 
ingly, the lower limits on have been dropped from the turn algorithm and the data 
given for case 28a of table V are without lower limits on kj. 

The work on the gain k^ showed that the critical value is the one determined at 
point A of figure 11 by equation (BIO). A correct value of k^ at this point insures a cor- 
rect value at point B and ensues a good landing. 


The Calculation of Gain kg 


The gain kg 



is calculated from the equation 
e 


(B12) 


and effectively gives the radians of bank per radian of error. The expression for kg 
shows that the only way in which kg can be changed is to vary k^. 

The effect of varying k^ on the touchdown conditions is presented in table IV and 
discussed in the subsection entitled "The Effect of Varying the Gain kg." A sizable 
group of initial conditions were run and an analysis of the results showed that the proper 
value of k^ could be calculated from 


k^ = 4.445 - 0.293;//. 


(B13) 


which makes the formulas for k^ 


kr 


/4.445 - 0.293;//A -j;//c-;// 


kl 


(B14) 


Remarks on the Turn Algorithm 

During the course of this investigation, several changes have been made to the turn 
algorithm as the result of wind disturbances and improved methods of determining the 
gains k| and kg. These changes are summarized in the turn algorithm shown in fig- 
ure 12. This turn algorithm is the one recommended for use with the autoland system 
described in reference 1. The improvement obtained through the use of this new algorithm 
is shown by a comparison of case 23 of table in (k^ = 4.4) with case 25a of table IV. 
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TABLE I.- CHARACTERISTICS OF TURBULENCE USED IN STUDY 



rms gust speed, 
m/sec, 
along axis - 

Maximum gust speed, 
m/sec, 
along axis - 

^1 

'^2 

^3 

^1 

^2 

^3 

— 
Mild turbulence 

0.163 

0.173 

0.185 

±0.45 

±0.45 

±0.55 

Medium turbulence 

1.59 

1.83 

1.88 

±4.5 

±4.6 

±5.5 

Heavy turbulence 

4.85 

5.5 

5.6 

±9.3 

±18.3 

±16.4 

Patchy turbulence 

1.177 

1.476 

1.353 

±5.86 

±7.23 

±9.73 



TABLE n.- SUMMARY OF THE EFFECT OF WIND, TURBULENCE, SAMPLE RATE, AND CONTROL-ACTUATOR NATURAL FREQUENCY ON THE AUTOLAND SYSTEM 


Tlnitial condition: = -6000 na; 


-4000 m; H = 540 m: >]/. = ' 


xr - 19 1A9 noa 


IS9 

CO 


Case 


Wind conditions 



System condition 

Measuring 

[ Airplane states 


1 





Average 

I 











no. 


Direction, 

rad 

Turbulence 

Shear 

Rotation 

Sample 

rate 





^2’ 

m 


«i. 


9, 

rad 


. 

Time, 

sec 



m/sec 

N.F., 

Hz 


I rad 

rad 

...... .. 1 


rad 

m/sec 

m/sec 

1 

0.0 

1 0.0 

— 

0.0 

None 

None 

1000 

30 

Start of flare 

3.82 X 10-3 

3.40 X 10-5 

-1.11' 

-3.14 

X ic2 

-5.75 X 10-2 

-8.16 X 10-2 

4.16 


116 

72.38 









Touchdown 

-9.41 X 10-4 

1.93 X 10-3 

-9.39 X 10-2 

2.84 

X 1q2 

-9.59 X 10-3 

-1.46 X 10-2 

6.71 X 

10-1 

124 

69.87 

2 ^ 

.0 

.0 

Light 

None 

None 

1000 

30 

Start of flare 

3.58 X 10-3 

3.23 X 10-4 

-1.06 

-2.98 

X 102 

-5.90 X 10-2 

-8.58 X 10-2 

4.30 


116 

72.93 









Touchdown 

-2.84 X 10-3 

8.76 X 10-4 

3.90 X 10"1 

2.65 

X 102 

-1.04 X 10-2 

-1.58 X 10-2 

7.32 X 

10-1 

124 

70.02 

3 b 

.0 

.0 

Medium 

None 

None 

1000 

30 

Start of flare 

3.03 X 10-3 

2.13 X 10-3 

-1.20 

-3.202 X 1Q2 

-7.09 X 10-2 

-1.05 X 10-1 

5.09 


115 

71.84 









Touchdown 

-1.3 X 10-2 

-9.62 X 10-4 

5.52 

2.31 

X lol 

-1.83 X 10-2 

9.98 X 10-3 

1.29 


120 

70.93 

4 b 

.0 

.0 

Heavy 

None 

None 

1000 

30 

Start of flare 

7.80 X 10-3 

2.00 X 10-4 

-1.28 

-3.31 

X 102 

-9.49 X 10-2 

-1.84 X 10-1 

1.52 


115 

83.21 









Touchdown 

2.43 X 10-2 

-1.80 X 10-2 

7.82 

7.39 

X lOl 

-5.68 X 10-2 

-1.14 X 10-1 

4.79 


118 

84.36 

5 

.0 

.0 

Medium 

None 

None 

1000 

30 

Start of flare 

5.94 X 10-3 

-7.12 X 10-3 

-1.22 

-2.49 

X 1q2 

-7.64 X 10-2 

-9.53 X 10-2 

5.52 


116 

72.31 









Touchdown 

3.74 X 10-2 

-8.26 X 10-3 

-1.25 

1.13 

X 1q2 

-6.29 X 10-3 

-6.48 X 10-3 

4.46 X 

10-1 

122 

71.08 

6 

.0 

.0 

Patchy 

None 

None 

1000 

30 

Start of flare 

1.69 X 10-2 

-1.09 X 10-2 

-2.80 

-3.32 

X 10^ 

-4.52 X 10-2 

-9.95 X 10-2 

3.28 


114 

72.81 









Touchdown 

-1.45 X 10-3 

-1.45 X 10-2 

-1.68 

9.69 

X lOl 

-6.72 X 10-3 

6.28 X 10-2 

4.94 X 

10-1 

120 

73.41 

7 ^ 

-25.80 

rr/4 

None 

None 

None 

1000 

30 

Start of flare 

9.98 X 10-3 

3.13 X 10-1 

-1.15 X to2 

-1.91 

X 102 

-2.62 X 10-2 

-1.06 X 10-1 

1.89 


111 

72.05 









Touchdown 

3.07 X lO-l 

5.56 X 10-1 

1.20 X 10-3 

1.32 

X 1q3 

-1.15 X 10-2 

-1.02 X 10-1 

7.79 X 

10-1 

133 

67.55 

a ^ 

-25.80 


None 

None 

None 

1000 

30 

Start of flare 

-1.69 X 10-2 

8.62 X 10-3 

2.63 

-2.77 

X 102 

-5.44 X 10--2 

-1.56 X 10-1 

3.98 


121 

73.11 









Touchdown 

-9.65 X 10-3 

-9.64 X 10-3 

2.16 

3.37 

X 102 

-8.18 X 10-3 

-1.15 X 10-1 

5.81 X 

10-1 

130 

71.06 

9 c,e 

-25,80 


None 

None 

None 

1000 

30 

Start of flare 

-5.73 X 10-3 

1.99 X 10-1 

-4.21 X lO-l 

-2.92 

X 10^ 

-6.06 X 10-2 

-1.51 X lO-l 

4.45 


103 

73.55 









Touchdown 

3.07 X 10-3 

1.97 X 10-1 

-3.08 X 10-1 

2.59 

X 102 

-9.74 X 10-3 

-1.04 X 10-1 

7.03 X 

10-1 

110 

72.15 

10 = 

-25.80 

7t/2 

None 

None 

None 

1000 

30 

Start of flare 

1.02 X 10-2 

3.43 X lO-l 

-5.35 

-3.10 

X 102 

-5.53 X 10-2 

-1.04 X 10-1 

4.07 


106 

73.64 









Touchdown 

5.91 X 10-3 

3.50 X 10-1 

-1.42 

3.23 

X lo2 

-8.11 X 10-3 

-4.60 X 10-3 

5.78 X 

10-1 

ll5 

71.35 

11 

-14.007 

7t/2 

None 

None 

None 

1000 

30 

Start of flare 

9.32 X 10-3 

1.93 X 10-1 

-3.52 

-2.97 

X 102 

-5.62 X 10-2 

-8.71 X 10-2 

4.07 


111 

72.54 







Touchdovm 

5.35 X 10-3 

2.01 X lO-l 

-1.36 

3.09 

X I02 

-8.54 X 10-3 

-2.37 X 10-2 

5.98 X 

10-1 

120 

70.05 

12 

-25.80 

i/4 

None 

Yes 

None 

1000 

30 

Start of flare 

-1.35 X 10-2 

2.43 X 10-2 

-1.22X 10-2 

-2.88 

X 102 

-6.03 X 10-2 

-1.19 X 10-2 

4.55 


100 

75.46 









Touchdown 

-7.32 X 10-3 

9.99 X 10-3 

-2.26 

9.67 

X lol 

-8,87 X 10-3 

4.61 X 10-2 

6.59 X 

10-1 

105 

74.26 

13 c.tl.e.g 

-25.80 

_ 

None 

None 

Yes 

1000 

30 

Start of flare 

4.85 X 10-2 

4.85 X 10-2 

-8.84 

-2.83 

X 102 

-5.84 X 10-2 

-1.85 X 10-2 

4.18 


111 

71.65 









Touchdown 

5.57 X 10-2 

1.14 X 10-1 

-7.22 

1.33 

X 102 

-1.38 X 10-2 

-1.79 X 10-2 

1.01 


117 

73.28 

14 <=A!.e 

-25.80 



None 

Yes 

Yes 

1000 

30 

Start of flare 

1.12 X 10-2 

1.30 X 10-3 

-2.59 

-2.83 

X 1q2 

-5.88 X 10-2 

-1.08 X 10-3 

4.31 


no 

73.40 









Touchdown 

8.48 X 10*4 

7.37 X 10-3 

-1.28 

1.51 

X 102 

-7.93 X 10-3 

6.17 X 10-2 

5.80 X 

10-1 

116 

73.18 

15 

-25.80 


Medium 

None 

None 

1000 

30 

Start of flare 

-9.69 X 10-3 

1.95 X 10-1 

3.23 

-3.09 

X 102 

-4.40 X 10-2 

-1.51 X 10-3 

3.25 


100 

73.85 









Touchdown 

3.46 X 10-3 

1.84 X 10-1 

2.04 

2.29 

X 102 

-1.13 X 10-2 

2.51 X 10-3 

8.50 X 

10-1 

107 

75.54 

16 

-25,80 


Medium 

Yes 

Yes 

1000 

30 

Start of flare 

-1.72 X 10-2 

-3.23 X 10-3 

6.10 

-3.05 

X 102 

-4.47 X 10-2 

-8.01 x'lO-2 

3.29 


100 

73.69 









Touchdown 

7.39 X 10-3 

-1.22 X 10-2 

-1.71 

4.80 

X 102 

-1.19 X 10-2 

1.18 X lO-l 

8.64 X 

lO-I 

111 

72.65 

17 c,d,e 

.0 

.0 

Medium 

None 

None 

100 

30 

Start of flare 

-6.52 X 10-3 

-9.35 X 10-3 

-2.63 

-3.09 

X 102 

-6.00 X 10-2 

-8.85 X 10-2 

4.28 


116 

71.36 









Touchdown 

2.26 X 10-2 

4.20 X 10-3 

-1.67 

1.94 

X 1q2 

-1.24 X 10-2 

5.37 X 10-2 

8.86 X 

10-1 

123 

71.36 

18 =■'> 

.0 

.0 

Medium 

None 

None 

50 

30 

Start of flare 

3.53 X 10-2 

1.28 X 10-2 

-1.55 X lOl 

-2.64 

X 1q2 

-7.07 X 10-2 

-1.22 X 10-1 

5.40 


104 

76.49 









Touchdown 

2.12 X 10-2 

2.54 X 10-2 

-9.26 

3.06 

X IQl 

-1.32 X 10-2 

3.43 X 10-2 

1.00 


108 

75.69 

19 

.0 

.0 

Medium 

None 

None 

10 

30 

Start of flare 

3.08 X 10-2 

1.21 X 10-2 

-1.38 X lOl 

-3.19 

X 1q2 

-7.35 X 10-2 

-1.12 X 10-1 

5.45 


104 

74.27 









Touchdown 

7.24 X 10-3 

9.80 X 10-3 

-8.04 

-1.67 

X lOl 

-2.84 X 10-2 

9.86 X 10-3 

2.09 


108 

73.58 

20 

.0 

.0 

Medium 

None 

None 

5 

30 

Start of flare 

5.81 X 10-3 

-7.94 X 10-3 

-1.46 

-2.74 

X lo2 

-7.27 X 10-2 

-1.29 X 10-1 

5.36 


116 

73.76 









Touchdown 

-1.19 X 10-2 

6.54 X 10-3 

-5.61 

3.67 


-3.19 X 10‘2 

3.92 X 10-2 

2.30 


119 

73.91 

21 c.d.e 

-25.80 

jt/4 

Patchy 

None 

None 

10 

30 

Start of flare 

-1.53 X 10-2 

1.98 X lO-l 

7.69 X 10-1 

-2.82 

X lo2 

-7.13 X 10-2 

-1.74 X 10-1 

5.19 


111 

72.92 









Touchdown 

-1.24 X 10-4 

1.86 X 10-1 

1.21 

2.13 

X I02 

-2.37 X 10-3 

-4.63 X 10’3 

1.71 


115 

72.15 

22 

-25.80 

ff/4 

Medium 

None 

None 

5 

5 

Start of flare 

5.10 X 10-2 

1.98 X 10-1 

-2.16 X iQl 

-2.76 

X 1Q2 

-4.12 X 10-2 

-1.44 X 10-2 

3.28 


99 

79.59 









Touchdown 

4.26 X 10-3 

2.19 X 10-1 

-9.85 

1.36 

X 102 

-1.56 X 10-2 

-2.77 X 10-2 

1.22 


104 

78.38 




















1 


= 77.12 for cases 1, 2, 3, 4, 5, 6, 17, 20, and 21; Vj = 78.56 for cases 13 and 14; Vj = 81.32 for case 8; Vj = 91.92 for case 11; = 97.09 for cases 7, 9, 12, 15, 16, 18, 19, and 22; Vj= 102.92 for case 10. 

b Basic system from NASA TN D-7611 (ref. 1). 

c This run has d feed to roll-control system and velocity compensation by equation (1), given by equation (2). 
d This run has special flap control for pitch angle. 

e The heading angle of the airplane in this run represents the crab angle with respect to the runway centerline. The velocity vector is pointing down the runway. 
fWind shear 0.0478 meter per second per meter. 

S Wind changed direction from tt/ 4 at 540 m to -?r/4 at touchdown. 




TABLE m.- THE EFFECT OF THE GAIN kj ON THE GROUND TRACK AND TOUCHDOWN CONDITIONS 
Initial condition; = -3000 m; ?2 = 4000 m; H = 540 m; i//. = 0.0 



Case 23 



Case 24 

Case 

* 

kj at A 

2.00 


1.796 


1.759 

kj atB* 

1.98 


1.753 


1.714 

T atB, sec 

28.6 


31.4 


31.7 

(p at B, rad 

1.039 x 10-3 


1.27 X 10-3 


9.60 X 10-4 

at B, rad 

1.29 


1,38 


1.39 

at B, m 

-1813.6 


-1825.1 


-1834.4 

?2 at B, m 

-2548.6 


-2327.5 


-2253.4 


Start of flare 

Touchdown 

Start of flare 

Touchdown 

Start of flare 

cp, rad 

-1.75 X 10-1 

-2.65 X 10-1 

-2.18 X 10-1 

-1.12 X 10-1 

-1.78X 10-1 

rad 

4.15 X 10-1 

2.0 X 10-^ 

2.97 X 10-1 

2.08 X 10-1 

2.61 X 10-1 

?2> “ 

-1.14 X 102 

-67.8 

-41.81 

-9.09 

-29.9 

m 

-2.48 X 102 

-73.3 

-334.1 

-40.9 

-334.4 

y, rad 

-1.43 X 10-1 

-6.86 X 10-2 

-9.82 X 10-2 

-2.61 X 10-2 

-7.43 X 10-2 

0, rad 

-2.46 X 10"^ 

-7.62 X 10-2 

-1.78 X 10-1 

-6.88 X 10-2 

-1.41 X 10-1 

W 2 j, m/sec 

10.45 

4.91 

7.85 

2.069 

5.41 

Time, sec 

71.0 

73.5 

73 

76.6 

74 

V, m/sec 

73.4 

7L6 



80.1 

79.3 

72.9 


-4.24 

158.8 


3.38 X 10-1 
80.9 
70.0 


Case 26 

1.706 I 

1.665 I 

32.0 j 

-8.87 X 10-4 j 
1.42 I 

-1851.0 ; 

-2252.2 I 

Start of flare Touchdown 

-1.24 x 10-1 -3.35 X 10-3 

2.21x10-1 1.90x10-1 

13.29 4.83 x 10-2 

-331.1 807.9 

-5.32 X 10-2 _5 30 x 10-3 

-1.18x10-1 9.95x10-2 

3.83 3.86 x 10-1 


See figure 11 for definition of these points. 


TABLE IV.- THE EFFECT OF VARYING THE GAIN k, ON TOUCHDOWN CONDITIONS, k, = ^ e 
[initial condition: Cj = 3000 m; Cg “ “4000 n\; H = 540 m.', 4 / = 0.0; as in cases 23 and 25 of table inj 



1 

Case 23 

Case 23a 

Case 23b 

Case 23c 

f 

Case 25 

Case 25a 

Case 25b 

‘'4 

4.8 


3.8 


3.2 


3.0 


4.8 


4.4 


4.0 



Start of flare 

Touchdown 

Start of flare ‘ 

Touchdown 

Start of flare 

Touchdown 

Start of flare 

Touchdown 

Start of flare 

Touchdown 

Start of flare 

Touchdown 

Start of flare 

Touchdown 

, (p, rad 

-1.75 X 10-' 

X 

o 

-2.28 X 10-' 

-2.78 X 10-' 

-2.70 X 10-' 

-3.11 X 10-' 

-2.88 X 10*' 

-2.79 X 10-1 

-1.78 X 10-1 i 

-1.95 X 10-2 

-1.89 X 10-1 

2.48 X 10-3 

-2,20 X 10-' 

3.32 X 10-2 

1 ^, rad 

4.15 X 10-' 

2.0 X 10-' 

3.53 X 10-1 

2.57 X to-' 

3.42 X 10-' 

2.29 X 10-1 

3.42 X lO-l 

2.11 X 10-1 

2.61 X 10-1 

2.10 X 10-1 

2.37 X 10-1 

1.96 X 10-' 

2.19 X 10-' 

1.69 X 10-' 

?2. I" 

-1.14 X 102 

-67.8 

-58.6 

-34.2 

-13.16 

-9.77 

47.8 

29.1 

-29.9 

-4.24 

-10.35 

1.94 X 10-' 

13.0 

4.1 

Cj, m 

-2.48 X 102 

-73.3 

• 243.5 

-128.2 

-294.2 

-170.4 

-321.5 

-171.9 

-334.4 

158.8 

-337.9 

196.4 

-349.5 

258.7 

Y, rad 

-1.43 X 10-> 

-6.86 X 10-2 

-1.42 X 10-' ' 

-7.32 X 10-2 

-1.43 X 10- 1 

-8.22 X 10-2 

, -1.42x10-1 

-7.07 X 10-2 

-7.43 X 10-2 

-4.83 X 10-3 

-7.24 X 10-2 

-7.47 X 10-3 

-7.12 X 10-2 

-1.0 X 10-2 

0, rad 

-2.46 X 10-' 

-7.62 X 10-2 

-1.88 X 10-' i 

-9.03 X 10-2 

-2.19 X 10-' 

-9.97 X 10-2 

-2.21 X 10-1 

-9.42 X 10-2 

-1.41 X 10-1. 

2.07 X 10-2 

-1,40 X 10-1 

3.02 X 10-2 

-1.48 X 10-' 

3.76 X 10-2 

Wjj, m/sec 

10.45 

4.91 

10.49 

5.43 

10.63 

6. 13 

' 10.39 

5.06 

5.41 

3.38 X 10-1 

5.29 

5.21 X 10-' 

5.36 

7.17 X lO-I 

: Time, sec 

Tl.O 

73.5 

72.0 1 

73.6 

72.0 

73.6 

72.0 

74.1 

74.0 

80.9 

74.0 

81.6 

74.0 

82.3 

V, m/sec 

73.4 

71.6 

74.1 

L 

74.3 

74.6 

1 

74.7 

73.4 

71.6 

72.9 1 

a 

70.0 

73.1 

. 

69.7 

75.3 

71.7 



TABLE V.- SUMMARY OF RESULTS OBTAINED USING NEW METHODS TO DETERMINE kj . AND k2 
[windspeed = -25.80 m/sec at ir/4 rad; patchy turbulence; sample rate 10 per seconc^ 



Case 27a 

Case 27b 

Case 27c 

Case 28a 

Case 28b 

Case 

28c 

?1. m 

-6000 


1 

-6000 


-6000 


-3000 


-3000 


-3000 


?2, m 

-4000 


-4000 


-4000 


-4000 


-4000 


-4000 


H, m 

540 


540 


540 


540 


540 


540 


1 //J, rad 

.0 


7t/4 


n/2 


.0 


7T/4 


n/2 



Start of flare 

Touchdown 

start of flare 

Touchdown 

Start of flare 

Touchdown 

Start of flare 

Touchdown 

start of flare 

Touchdown 

Start of flare 

Touchdown 

rad 

-5.64 X 10-3 

2.36 X 10-3 

-1.10 X 10-3 

2.49 X 10-2 

3.04 X 10-3 

-3.66 X 10-3 

-1.51 X 10-1 

1.52 X 10-2 

-7.89 X 10-3 

1.91 X 10-3 

-8.00 X 10-3 

-1.04 X 10-3 

x^, rad 

1.95 X 10-1 

1.86 X 10-1 

1.95 X 10-1 

1.97 X 10-1 

1.94 X 10-1 

1.85 X 10-1 

2.03 X 10-1 

1.87 X 10-1 

2.04 X 10-1 

1.99 X 10-1 

2.32 X 10-1 

1.84 X 10-1 

Cj, m 

4.25 X 10-1 

4.18 X 10-1 

1.23 

-1.27 X 10-1 

-6.46 X 10-3 

6.82 X 10-1 

1.51 

3.66 X 10-1 

-4.93 X 10-1 

-2.88 X 10-1 

1.57 X lol 

-1.05 

Cj, m 

-2.91 X 10^ 

2.32 X 103 

-2.88 X 1q3 

1.37 X lo3 

-2.96 X lo3 

8.45 

-3.02 X lo3 

3.44 X 103 

-3.36 X lo3 

5.18 X IqI 

-3.51 X 103 

8.72 X 1q3 

y, rad 

-4.12 X 10-3 

-9.93 X 10-3 

-5.65 X 10-3 

-1.63 X 10-3 

-6.23 X 10-3 

-7.79 X 10-3 

-4.74 X 10-3 

-1.15 X 10-3 

-7.70 X 10-3 

-1.41 X 10-3 

-6.31 X 10-3 

-2.61 X. 10-3 

9, rad 

-1.5» X 10-1 

1.60 X 10-2 

-1.70 X 10-1 

-9.56 X 10-3 

-1.70 X 10-1 

-3.33 X 10-3 

-1.19 X 10-1 

5.68 X 10-3 

-1.40 X 10-1 

5.86 X 10-3 

-1.38 X 10-1 

1.04 X 10-1 

W 2 j, m/sec 

2.98 

7.27 X 10-1 

4.06 

1.18 

4.50 

5.68 X 10-1 

3.54 

8.34 X 10-1 

5.61 

9.83 X 10-1 

4.78 

1.95 X 10-1 

Time, sec 

114.0 

121.7 

117.0 

122.8 

116.0 

120.2 

75.0 

83.8 

74.0 

79.5 

75.0 

91.3 

V, m/sec 

72.28 

73.24 

71.91 

72.45 

72.24 

72.94 

74.69 

72.49 

72.97 

69.65 

75.86 

74.7 



Enhanced 

reflector 

at 

runway 


Ground 

radar 


Notes: 


Airborne radar data 


Line of sight 


Sensed airplane motions 


Airborne 

radar 


W/Z/Z/A, 


Line of siaht 


Airborne 

processor 


Wind and turbulence 




(1) 77/7/97 Indicates equipment fixed to airplane 

that senses airplane motions 

(2) IMU, Inertial measuring unit 

(3) a,fl,D,and & are computed in the IMli 

<4)aii actuators are rate and displacement limited 


SnT 


'Tc- 


[Coordinate! 

trans- 

jformation 




Ground radar data 




Turn 

computer 




Let down 
comouter 


Input, output of this 
section varied for 
sample rale 


‘t’c 


.6^= kgP- k,p- 6^- C^(3- kjjE 


Roil autopilo* / ■ 

6 ^ . 

\ MVV 


l+rD 

T=0. 03779, 


^0- P.P 


Let down autopilot | 6g 


Servo 

1 

1+tD 

r*0, 039641 


Spoilers 


Stabilizer autopilot 
1 Cm^P 
5h‘ 


Servo 
1 

1 +tD 

T*0. 03779 


P.P 


Speed control 
2^7' (V^-Vi-mu 
,k5=Vc(28.13I+58, 9736fi 


Thrust change 


-Forces and moments 


Airplane motions 


Forces 

and 

moments } 


^0 

Rudder autopilot 


[ Servo 
! 1 

u 

1 ‘^0 T 


6, ■ -kjjp-ki3& 

1 

'l +rD 

1 



1 

1 ^ -0. 03964 



Figure 1.- Block diagram of the automatic landing system. Symbols are defined in reference 1. 
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Speed, Wyy , m/ sec Speed, Wyy , ml sec 



Time.t, sec 



Time, t, sec 
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(c) Turbulence along ?72-axis 
Figure 2 ,- Concluded. 



Roll angle, f, rad Roll angle, <p , rad Roll angle, <p, rad 











Pitch angle, 0 , rad Pitch angle, e,rad Pitch angle, 0 ,rad 


2 




Time, sec 


(c) Pitch angle. 
Figure 3,- Continued. 
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no turbulence, case 1 


I I ! I M I I I I I i I I ! ! I : I ' M I j i I ! ; I i J I ! I ! I I I i i 

20 40 60 80 

Time, sec 


medium turbulence, case 3 " 





medium turbulence, case 3 


} 40 60 80 

Time, sec 


patchy turbulence, case 6 . 

A , 

' V 



Time, sec 

(e) Flight -path angle, 
Figure 3.- Concluded 






Altitude, H, m 


Co 

00 



(a) Altitude track. 

Figure 6.- Ground and altitude tracks for the initial condition = -6000 m; = -4000 m; H = 540 m; 

= 7t/2; windspeed -25.80 m/sec at tt/ 4; and medium turbulence. 




S|,m 

(b) Ground track. 
Figure 6.- Concluded. 


CO 

CO 




Figure 8.- Comparison of groiind tracks for correct, case 25, and incorrect, case 23 

settings of k^. 






Figure 9. - Comparison of ground tracks for final method of determining kj and kg. Patchy turbulence and wind 
-25.80 m/sec at 7r/4 rad. Details of initial condition and final conditions are given in table V. 



H3 


Notes 


(1) Runway coordinates are centered at touchdown point ; positive directions 
shown at right 

(2) Airplane inertial coordinates are/' to runway coordinates but centered 
in airplane 


(3) The Z| coordinates are centered in the airplane and rotated away fro 

the airplane inertial coordinates by the Euler angle i/i the heading angle 
of the airplane. The following sketch shows these relationships; 



(4) For the runway i/i = i/i^ = 0 

(5) The '1\'V2 plane coincides 
with the Z|, z^ plane 


Airplane 

location 





Figure 10.- Geometry of a typical landing situation 


that shows coordinate systems used in study. 



A First determination of k 



Figure 11,- Sketch of ground track when = 0. 




Notes:- (1) This diagram assumes zero roll angle at the start of the landing maneuver 

(2) BLIP, CLIP, and CHIP are switch variables and are either true (T) or false (F) 


Figure 12.- Turn-control algorithm. Equation numbers refer to equations in the text. 
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